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Abstract 

We  carry  out  an  energy-inventory  analysis  for  a  12-cm- 
diameter,  argon  gas-puff  shot  on  the  Saturn  generator. 
From  the  measured  pinch  inductance  time  history,  we 
infer  the  average  radius  of  the  pinch  current,  rind,  and  the 
total  work  done  on  the  pinch,  Ecoupied.  both  as  functions  of 
time.  The  kinetic  energy  as  a  function  of  time.  Eke,  is 
inferred  from  rind  and  the  measured  initial  mass 
distribution,  assuming  all  of  the  mass  is  accreted  as  in  a 
snowplow.  The  results  show  that  early  in  the  implosion, 

Ecoupied  Eke-  However,  when  the  pinch  begins  to  radiate 
significantly.  Eke  begins  to  decrease  while  Ecoupied 
continues  to  increase.  The  increase  in  coupled  energy 
goes  into  plasma  internal  energy.  Hint,  i.e.,  ionization  and 
thermal  energies  of  the  z-pinch  plasma,  and  total 
radiation,  Ej-ad-  During  the  time  of  K- shell  emission,  the 
significant  (~  150-kJ)  increase  in  Ecoupied  is  reflected  in 
approximately  the  same  increase  in  Ej-ad-  There  is 
agreement  between  Hint  at  the  time  of  peak  K-shell 
emission  inferred  from  Ecoupied  and  measurements  of  Ej-ad, 
and  Hint  inferred  from  independent  spectroscopic 
measurements  combined  with  an  atomic-physics  model. 

I.  Introduction 

We  wish  to  maximize  the  K-shell  yield  from  high- 
atomic-number,  gas-puff  z  pinches,  e.g.,  K-shell  energy 
of  ~  1  -  8  keV,  and  to  develop  z-pinch-based  concepts  for 
sources  of  higher  photon  energy  (up  to  40  keV).  To  this 
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end,  large-diameter  implosions  [1-8}  are  needed  to:  (a) 
produce  the  high  specific  energy  required  to  excite  K- 
shell  line  radiation  from  high-atomic-number  z  pinches, 
(b)  properly  match  the  load  to  high-current  (~  10  MA  and 
higher)  generators  [9-12],  and  (c)  develop  continuum- 
radiator  concepts  [13-15]  that  might  require  ionizing 
beyond  the  optimum  He/H-like  state  for  K-shell  emission. 
Large-initial-diameter  implosions  are  also  required  for 
efficient  coupling  between  the  z  pinch  and  a  generator 
with  a  long  current  rise  time  (>  100  ns)  and,  thus,  a 
desirable  lower  driving  voltage.  It  has  been  shown  for 
gas-puff  z-pinches  both  theoretically  [16-18]  and 
experimentally  [7,19-21]  that  the  expected  deleterious 
effects  on  final  pinch  formation  associated  with  the 
increase  in  the  Rayleigh-Taylor  instability  growth  at 
larger  diameter  can  be  mitigated  by  using  an  initial  mass 
profile  that  is  peaked  on  axis. 

An  important  aspect  of  z-pinch  physics  is  obtaining 
an  understanding  of  the  pinch  energetics.  For  example, 
work  with  wire  arrays  suggested  that  more  energy  is 
radiated  than  can  be  accounted  for  by  one-dimensional 
implosion  calculations  of  kinetic  energy  that  assume  the 
mass  resides  in  a  thin  shell.  [22]  A  higher-than-expected 
resistivity  was  one  proposed  explanation  for  this 
observation.  As  noted  in  previous  simulations  of  wire 
arrays,  the  total  work  done  on  the  pinch  plasma  by  the 
Lorenz  force  can  far  exceed  the  kinetic  energy. [23] 
Investigations  into  the  energetics  of  wire  arrays 
continue.  [24,25] 
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For  the  gas-puff  implosions  with  structured  radial 
profiles  studied  here,  the  thin-shell  approximation  is 
clearly  not  appropriate.  Moreover,  the  fact  that  the  initial 
mass  is  distributed  in  radius  also  calls  into  question  the 
role  of  kinetic  energy  in  the  overall  pinch  energetics,  as 
much  of  the  mass  is  not  initially  at  large  radius.  The 
transfer  of  energy  from  ions  to  electrons  during  the 
stagnation  of  a  neon  gas-puff  z-pinch  has  been 
studied.  [26]  However,  the  global  energy  inventory  for 
large  diameter  gas-puff  z  pinches  has  not  been  explored. 

In  this  paper,  we  focus  on  Saturn  shot  3565,  which 
was  one  of  several  taken  during  the  first  experiments  on 
Saturn  with  a  12-cm  diameter,  triple-shell,  argon  gas  puff. 
This  particular  shot  was  chosen  because  all  the  diagnostic 
signals  were  available  and  a  full  analysis  was  possible. 
Also,  this  was  the  only  shot  in  the  series  for  which  the 
vacuum  voltmeter  was  well  calibrated  and  properly 
shielded,  [27]  allowing  for  the  first  time  a  reliable 
determination  of  the  voltage  driving  the  gas  puff  as  a 
function  of  time.  As  will  be  seen,  this  is  a  critical 
measurement  for  determining  the  overall  energy 
inventory.  From  this  inventory,  the  plasma  internal 
energy,  which  is  the  sum  of  the  thermal  and  ionization 
energies  of  the  pinch  plasma,  is  inferred.  The  internal 
energy  so  inferred  is  compared  to  the  internal  energy 
obtained  independently  from  spectroscopic  measurements 
and  an  atomic  physics  model.  Knowledge  of  the  initial 
mass  distribution  also  is  critical  for  these  analyses. 

The  details  of  the  experimental  configuration  and 
description  of  the  energy-inventory  analysis  are  given  in 
Section  II.  In  Section  III  we  present  data  and  analyses  for 
the  overall  energy  inventory  of  Saturn  shot  3565.  In 
Section  IV,  the  internal  energy  near  the  time  of  peak  K- 
shell  power  obtained  from  the  energy  inventory  is 
compared  to  that  obtained  from  the  measured  spectral 
characteristics  in  combination  with  an  atomic  physics 
model.  The  work  is  summarized  in  Section  V. 


II.  EXPERIMENTAL  SETUP  AND 
DESCRIPTION  OF  ANALYSIS 

The  work  reported  herein  was  performed  on  the 
Saturn  generator, located  at  Sandia  National 
Laboratories.  For  this  work,  the  posthole  convolute 
diameter  was  increased  from  15  cm  to  30  cm.^^  Also,  the 
Saturn  pulse-forming  line  was  shorted  out  to  provide  the 
“long-pulse”  mode,  where  the  current  rise  time  into  a 
short  circuit  load  is  ~  230  ns.  The  gas-puff  nozzle  is  a  12- 
cm  diameter,  triple- shell  nozzle  illustrated  in  Fig. 
I. [7, 19,20]  The  initial  gas  mass  distribution  is  measured 
using  the  planar,  laser-induced  fluorescence  (PLIF) 
method. [30,3 1]  For  Saturn  shot  3565,  the  three  plenums 
were  pressurized  as  follows:  outer  =  3.3  psia,  inner  =  lO.O 
psia,  center  jet  =  72.0  psia.  This  combination  resulted  in 
nominal  average  mass/length  of  425  pg/cm  with  an 
uncertainty  of  ±15%  and  a  radial  distribution  that  is 


peaked  on  axis. [7, 19]  The  latter  is  essential  for  achieving 
high-quality  implosions  from  large  radius.  The  radiation 
diagnostics  include  a  bare  bolometer  (nickel  element, 
corrected  for  K-transparency)  to  measure  the  total 
radiated  energy,  filtered  calorimeters  to  measure  the  K- 
shell  yield,  and  a  variety  of  photo-conductive  devices 
(PCDs)  filtered  to  measure  the  argon  K-shell  time  history. 
For  this  work,  the  PCDs  were  normalized  to  the 
calorimeters.  A  K-shell,  time-integrated  image  of  the 
pinch  is  obtained  from  an  appropriately  filtered  CCD- 
based  pinhole  camera.  Argon  K- spectra  and  the  K- 
spectra  from  a  chlorine  dopant  (Freon),  introduced  in 
either  the  center  jet  or  inner  shell  on  various  shots,  are 
measured  using  a  CCD-based  crystal  spectrograph.  The 
L- shell  radiation  is  inferred  from  a  filtered  x-ray  diode 
normalized  to  a  bolometer. 

The  current  is  determined  by  subtracting  and 
averaging  the  signals  from  two,  oppositely  wound  B-dots 
180°  apart  in  azimuth  and  located  between  the  posthole 
convolute  and  the  gas-puff  z-pinch.  The  voltage  driving 
the  gas  puff  is  directly  measured  using  a  vacuum 
voltmeter[32]  The  implementation  of  this  diagnostic  on 
Saturn  is  described  in  detail  in  Ref  27. 


center  jet 


Figure  1.  Cross  section  of  the  12-cm 
diameter,  triple  nozzle  used  on  Saturn 
shot  3565. 


Figure  2  gives  the  current  and  K-shell  x-ray  power 
as  functions  of  time  for  Saturn  Shot  3565.  The  peak 
current  is  6.5  MA  and  the  implosion  time  is  206  ns.  The 
energy  radiated  from  the  K-shell  for  this  shot  is  72  kJ, 
roughly  twice  what  was  previously  obtained  at  shorter 
implosion  time  (~  50  ns)  using  a  smaller-diameter  (~  few 
cm)  gas-puff. [20,33] 

The  total  energy  delivered  to  the  z-pinch  plasma 
from  the  electrical  circuit,  Ecoupied.  is  given  by^^ 

(1). 
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where  I  and  V  are  the  z-pinch  current  and  driving  voltage, 
respectively,  and  L  is  the  z-pinch  inductance.  In  Eq.  (1), 
we  have  assumed  that  that  resistance  is  negligible.  The 
change  in  L,  AL,  can  be  expressed  as: 


f 

lSL{nH)  =  L{t)-L^  =  2x/(cm)xln 

V 


(2). 


Here  Lq  is  the  fixed  inductance  associated  with  the  region 
between  the  location  where  the  voltage  is  measured  and 
the  initial  z-pinch,  I  is  the  pinch  length  (3.8  cm  for  this 
work),  rind  is  the  axially- averaged  effective  radius  of  the 
pinch  current,  and,  again  assuming  resistance  is 
negligible,  L(t)  is  given  by 


i(r)= 


(3). 


Ecoupied  is  the  total  work  done  on  the  plasma  and  can  be 
expressed  as 


^coupled  (0  “  ^KE  ^int  ^ rad 


where  Eke  is  the  kinetic  energy  associated  with  the 
radially-directed  pinch  implosion  and  Ej-ad  is  the  total 
radiated  energy.  The  internal  energy,  Eint  is  given  by 

(5), 

where  Z  is  the  is  the  average  ionization  state  of  the  pinch 
plasma,  Ni  is  the  total  number  of  ions  in  the  pinch,  k  is 
Boltzman’s  constant,  Ti  (Te)  and  is  the  ion  (electron) 
temperature,  and  Eionize  (a  well  known  function  of  Z)  is 
the  energy  required  to  ionize  up  to  the  ionization  state  Z. 
The  second  term  on  the  right-hand  side  of  Eq.  (5) 
represents  the  plasma  thermal  energy  associated  with 
random  motion.  Eor  this  work,  EKE(t)  is  obtained  using 
d(rind)/dt  for  the  radial  speed  and  the  measured  initial 
mass  distribution  to  determine  the  mass  accreted  up  to 
rind(t)  with  d(rind)/dt,  so  that 


^KE  (0  “  2  ^accret  (0 


dr, 


ind 

dt 


(6), 


Figure  2  Plot  of  current  and  K-shell  power  as 
functions  of  time. 


where  maccret  is  the  total  mass  accreted  up  to  at  time  t. 


Ill  MEASUREMENTS  AND  ANALYSIS 
FOR  ENERGY  INVENTORY 

Plotted  in  Eigure  3  are  the  pinch  current,  K-shell 
power,  and  rind  obtained  from  Eq.  (2),  all  as  functions  of 
time.  Also  shown  are  two  radii  from  a  0-D  snowplow 
calculation,  rgp,  using  the  measured  current  and  measured 
initial  mass  distribution  with  two  total  masses,  361  pg/cm 
(within  the  measurement  uncertainty)  and  425  pg/cm 
(nominal).  In  general,  the  snowplow  calculation  suggests 
that  the  independently  measured  mass,  current,  and 
implosion  time  are  consistent.  Using  361  pg/cm  gives  a 
better  match  to  the  206-ns  implosion  time.  Also  note  that 
the  rsp  is  continually  accelerating  through  stagnation  while 
rind  decelerates.  This  is  probably  a  consequence  of  the 
snowplow  calculation  not  accounting  for  the  plasma 
pressure.  Unlike  rgp,  the  implosion  speed  given  by 
d(rind)/dt  is  remarkably  constant,  46  cm/ps.  This 
constant  speed  implies  that,  on  average,  the  radial 
acceleration  is  very  small  (a  result  of  the  initial  radially 
distributed  mass  distribution),  suggesting  that  the 
Raleigh-Taylor  growth  is  mitigated.  The  initial  radius  for 
the  snowplow  calculation,  8.3  cm,  was  chosen  based  on 
the  largest  radius  for  which  there  are  PLIE  data.  More 
sensitive  measurements  have  shown  the  presence  of  gas  at 
large  radii  in  similar  nozzles. [34]  The  8. 3 -cm  radius  is 


time  (s) 

Figure  3  Plot  of  pinch  current,  K-shell  power  and  rind 
from  Eq  (2).  Also  shown  are  two  radii,  rgp,  calculated 
using  the  0-D  snowplow  model  with  the  pinch  current 
and  measured  initial  density  distribution  for  two 
different  total  masses  as  inputs. 
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PK-shell  (TW) 


very  close  to  the  initial  value  of  rind-  However,  both  are 
very  close  to  the  return  current  post  radius  of  8.6  cm.  All 
of  this  suggests  that  the  initial  current  path  might  not  be 
purely  axial,  but  have  a  significant  radial  component  that 
results  in  the  initial  current- sheath  motion  having  a 
component  in  the  axial  direction. [35] 

Ecoupied  from  Eq.  (1)  along  with  Eke  from  Eq.  (5)  are 
shown  in  Eig.  4  along  with  the  pinch  current  and  K- shell 
power,  all  as  functions  of  time.  Early  in  the  implosion, 
Ecoupied  ~  Eke-  Because  the  implosion  speed  is  nearly 
constant,  the  increase  in  Eke  is  a  result  of  mass  accretion. 
At  rind  ~  2.3  cm,  when  Ecoupied  -150  kJ,  Eke  begins  to 
decrease  as  the  pinch  begins  to  slow  down.  This  is  about 
the  time  the  total  radiation  becomes  appreciable  (see  Eig. 
5).  The  rate  of  slowing  down  more  than  compensates  for 
the  small  but  finite  mass  accretion  at  small  radius.  As 
Eke  begins  to  decrease,  however,  Ecoupied  continues  to 
increase  over  the  next  -  70  ns  to  440  kJ  at  the  time  of 
peak  K- shell  power.  Also  plotted  in  Eig.  4  is  Ecoupied  - 
Eke  =  Eint  +  Erad  (see  Eq.  [4]).  Near  the  time  of  peak  K- 
shell  power.  Eke  has  decreased  nearly  to  zero  and  Ecoupied 

“  Ejnt  +  Ej-ad- 

The  total  radiated  energy,  E^d,  is  plotted  on  an 
expanded  time  scale  in  Eig.  5,  along  with  the  pinch 
current,  K- shell  power,  and  Ecoupied-  During  the  time  of 
K-shell  emission,  Ecoupied  increases  by  150  kJ,  from  350  kJ 
to  500  kJ,  showing  that  energy  continues  to  be  delivered 
to  the  pinch  during  stagnation.  During  that  same  time 
interval,  the  total  radiation  increases  by  nearly  the  same 
amount,  140  kJ,  and  the  total  K-shell  yield  is  72  kJ.  A 
significant  increase  in  coupled  energy  during  the  time  of 


time  (s) 

Figure  4  Plot  of  Ecoupied  and  Eke  as  functions  of 
time  along  with  pinch  current,  K-shell  power  and 
rind  all  as  functions  of  time.  At  time  of  peak  K-shell 
power,  Ecoupied  =  Eint  +  Erad 
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Figure  5  Plot  of  Ecoupied  and  Emd  as  functions  of 
time  along  with  pinch  current  and  K-shell  power. 
During  the  time  of  K-shell  emission,  Eint  ~  Ecoupied  - 
Erad  ~  constant,  250  kJ. 


K-shell  emission  also  was  observed  in  a  neon  gas-puff  z- 
pinch  pinch  where  a  vacuum  voltmeter  was  first  used  to 
obtain  the  coupled  energy. [36]  As  noted  above  (see  Eig. 
4),  Eke  drops  to  near  zero  during  the  K-shell  emission. 
Thus,  from  Eq.  (4),  Eint  -  Ecoupied  -  Erad-  and  as  seen  in  Eig. 
4,  Eint  is  constant  at  about  250  kJ,  during  this  time.  Thus 
a  steady  state  of  sorts  is  established  during  the  time  of  K- 
shell  emission,  with  the  coupled  energy  going  into 
radiated  energy  and  the  internal  energy  remaining  fixed. 
The  process  is  not  direct,  however.  It  is  thought  that  the 
higher  temperature  ions  transfer  energy  to  electrons  which 
then  lose  their  energy  via  radiation. 

In  the  next  section,  we  compare  Eint  determined  from 
electrical  and  total  radiation  measurements  at  the  time  of 
peak  K-shell  power  (250  kJ),  with  the  internal  energy  at 
that  time  obtained  from  independent  spectroscopic 
diagnostics  in  combination  with  an  atomic  physics  model. 


IV  SPECTROSCOPIC  DETERMINATION 
OF  INTERNAL  ENERGY 

Using  a  well  known  technique, [37,38]  the  argon  K- 
shell  spectra  can  be  used  to  obtain  the  radially-integrated 
axially-averaged  Te.  The  argon  atomic  model  used  in  this 
technique  has  been  improved  and  now  includes  detailed 
level  structure  for  all  1 8  ionization  stages  of  argon,  a  total 
of  185  atomic  levels.  The  radiation  transport  equations 
are  solved  for  611  lines  with  10  to  25  photon  energy 
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points  per  line.  The  same  procedure  is  applied  to  the 
chlorine  K-shell  spectra  to  obtain  a  radially  localized  Te 
near  the  time  of  peak  K-shell  power  in  a  central  region 
(Freon  in  the  center  jet)  or  an  inner  region  (Freon  in  the 
inner  shell).  Knowledge  of  Tg,  combined  with  the 
absolute  K-shell  power  and  measurement  of  the  K-shell 
emitting  region  is  used  to  obtain  the  ion  density  near  the 
time  of  peak  K-shell  power,  ni.[37,38]  This  analysis 
reveals  for  Saturn  shot  3565  that  Tg  in  the  central  and 
inner  regions  is  1.9  and  1.4  keV,  respectively.  In 
addition,  44%  of  the  initial  mass  is  radiating  in  the  K- 
shell  from  a  region  of  1.15  mm  in  radius,  with  Ui  =  7x10^^ 
cm'^. 

Based  on  these  measured  quantities,  we  construct  a 
3 -radial-zone  model  for  the  stagnated  pinch  plasma  in 
which  the  two  central  zones,  zone  1  (central)  and  zone  2 
(inner)  radiate  in  the  K-shell  with  two  different  Te’s  and 
with  mass  ratio  given  by  the  measured  initial  injected 
mass  ratio  between  the  nozzle  center  jet  and  the  inner 
shell,  1:2.  The  two  central  zones  are  surrounded  by  a 
colder  blanket,  zone  3,  that  comprises  the  remainder  of 
the  mass,  56%  of  the  initial  mass,  and  radiates  only  in  the 
L-shell.  The  measured  L-shell  power  at  the  time  of  peak 
K-shell  emission,  4  TW/cm,  is  used  with  the  argon  atomic 
physics  model  to  infer  the  ion  density  and  thus  the  radius 
associated  with  the  L-shell  radiating  region.  We  assume 
the  electron  temperature  in  the  outer  zone  is  0.7  keV  (half 
of  Tg  in  zone  2).  We  have  no  measurement  of  Ti,  so  we 
assume,  based  on  previous  z-pinch  experiments  where  Ti 
was  measured,  that  Ti  ^  20  keV  in  all  three  zones.  Note 
that  if  one  equates  the  ion  kinetic  energy  to  a  thermal 
energy  of  3Ti/2,  assuming  the  observed  46  cm/ps  radial 
velocity,  one  obtains  Ti  =  28  keV  as  an  upper  bound  to  Ti. 
Thus,  the  assumption  for  Ti  is  quite  reasonable.  The 


Figure  6  Schematic  and  table  for  the  three-zone 
model.  Parameters  in  orange  are  not  inferred  from 
measured  quantities. 


atomic  physics  model  self-consistently  determines  Z  for 
the  given  plasma  conditions. 

Figure  6  summarizes  the  three-zone  model, 
measurements,  and  assumptions  described  above. 
Recalling  Eq.  (5)  for  Eint,  we  can  see  from  Fig.  (6)  that  the 
total  internal  energy  can  be  calculated  from  summing 
across  the  three  zones.  This  method  gives  252  kJ  for  the 
total  internal  energy  at  the  time  of  peak  K-shell  emission 
compared  with  250  kJ  determined  from  the  energy 
inventory  (see  Fig.  5  and  accompanying  text).  While 
there  are  many  assumptions,  this  is  quite  remarkable 
agreement.  In  future  experiments,  attempts  will  be  made 
to  reduce  the  uncertainties  in  this  approach. 


III.  SUMMARY 

We  have  carried  out  an  energy  inventory  analysis 
for  a  12-cm-diameter,  argon  gas-puff  shot  on  the  Saturn 
generator  with  a  peak  current  of  6.5  MA  and  an  implosion 
time  of  206  ns.  This  resulted  in  72  kJ  of  argon  K-shell 
radiation,  about  a  factor  of  2  more  than  previously 
obtained  on  Saturn  using  smaller  diameter  nozzles  and  at 
shorter  implosion  times.  Using  the  measured  current  and 
voltage,  we  determine  Egoupied.  the  total  work  done  on  the 
pinch  plasma,  and  the  pinch  inductance  time  history,  from 
which  we  infer  the  axially  averaged  radius  of  the  pinch 
current,  qnd,  both  as  functions  of  time.  The  kinetic  energy 
as  a  function  of  time.  Eke,  is  inferred  from  rind  and  the 
measured  initial  mass  distribution,  assuming  all  the  mass 
is  accreted  as  in  a  snowplow.  The  results  show  that  early 
in  the  implosion,  Ecoupied  ~  Eke-  However,  when  the 
Ecoupied  ~  150  kJ,  rind  ~  2.3  cm,  and  the  pinch  begins  to 
radiate  significantly.  Eke  begins  to  decrease  while  Ecoupied 
continues  to  increase.  Over  the  next  ~  70  ns.  Eke 
becomes  negligible  while  Ecoupied  increases  to  440  kJ  at 
the  time  of  peak  power.  The  increase  in  coupled  energy 
goes  into  Eint  and  Emd-  During  the  time  of  K-shell 
emission,  Ecoupied  increases  by  150  kJ,  showing  that 
significant  amount  of  energy  continues  to  be  delivered  to 
the  pinch  during  stagnation.  This  increase  in  Ecoupied  is 
reflected  in  approximately  the  same  increase  in  Emd,  with 
Eint  ~  constant  (250  kJ),  suggesting  a  steady  state  where 
Ecoupied  is  rapidly  radiated  away.  There  is  agreement 
between  the  value  Eint  at  the  time  of  peak  K-shell 
emission  inferred  from  Ecoupied  and  measurements  of  Emd, 
and  Eint  inferred  from  independent  spectroscopic 
measurements  combined  with  an  atomic  physics  model. 
This  agreement  helps  to  confirm  our  approach  to  the 
energy-inventory  analysis.  It  is  interesting  to  note  that 
because  wire  arrays  rapidly  transform  from  a  line  mass  to 
a  distributed  mass,  the  conclusions  reached  in  this  gas- 
puff  work  may  apply  to  wire  arrays  as  well. 

Measurement  of  the  pinch  driving  voltage  is  an 
essential  part  of  carrying  out  this  energy  inventory.  We 
hope  to  apply  this  approach,  with  improved  current  and 
radiation  diagnostics,  to  future  experiments  on  Saturn. 
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